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ABSTRACT 

bo" 
on . 

Q\ , We have used spectra taken between 1992 and 1997 to derive the spectroscopic 

orbit of the eclipsing double-lined spectroscopic binary HD 197770. This binary has 
a period of 99.69 ± 0.02 days and K amplitudes of 31.2 ± 0.8 and 47.1 ± 0.4 km s" 1 

■ for components A & B, respectively. The msin 3 i values for A & B are 2.9 and 1.9, 
O ■ respectively, and are close to the actual masses due to the eclipsing nature of this 

binary. Both components of HD 197770 have spectral types near B2 III. This means 
both components are undermassive by about a factor of five and, thus, evolved stars. 
Additional evidence of the evolved nature of HD 197770 is found in 25, 60, and 
| 100 fim IRAS images of HD 197770. These images show 2 apparent shells centered 

on HD 197770; a bright 60 /jm shell with a 14' diameter and a larger (1?2 diameter) 

■ bubble-like feature. At least one of the components of HD 197770 is likely to be a 
post-AGB star. 
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Introduction 



Interest in the star HD 197770 (HR 7940, a(2000) = 20 h 43 m 13!52, 5(2000) = +57 o 6'50"9) 
increased greatly with the discovery that its line-of-sight has a polarization feature coincident with 



the 2175 A extinction bump ( piayton et al. 1992 ; Anderson et al. 1996). Out of the 30 sightlines 



with UV spectropolarimetry, such a polarization bump has only been seen along one other sightline 
( Wolff et al. 199^ ). HD 197770 has long been known to have a variable radial velocity (Adams 



Joy, fc Sanford 1924 ; piaskett fc Pearce 193CQ . Observat ions originally intended to study the 
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sightline towards HD 197770 have shown it to be a double- lined spectroscopic binary (Hanson 
|et al. 1994 ; piayton 1996| ). Recent photometric observations have shown that HD 197770 is an 
eclipsing binary ( Jerzykiewicz 1993| ; piayton 1996 ). We obtained spectra of HD 197770 between 
1992 and 1997 in order to determine the spectroscopic orbit of HD 197770. The spectroscopic 
orbit, coupled with the eclipsing nature of the binary, allowed us to determine the masses of the 
binary components. 



2. Observations 

In 1992 June and October, observations of HD 197770 were acquired at the Kitt Peak 
National Observatory (KPNO) 0.9m coude feed telescope using the coude spectrograph. The June 
observation was obtained using an echelle grating and a cross disperser grism to give disjoint orders 
between 3870 and 4085 A with a resolution of 120,000. It was reduced using IRAF. The October 
observations also used an echelle grating and a cross disperser grism giving orders which covered 
5580-7160 A (HJD 2448901.6) or 3990-4490 A (HJD 2448902.6 & 2448903.6) at a resolution of 
80,000. Between 1996 September and 1997 October, 17 observations of HD 197770 were taken 



at the lm Ritter Observatory telescope using a fiber-fed echelle spectrograph (Gordon Mulliss 



1997). These observations cover nine disjoint 70 A-orders between 5200 and 6600 A at a resolution 



of 25,000. The KPNO and Ritter observations were reduced using a package written to reduce 
Ritter observations. Since the KPNO observations were not acquired with a fiber-fed echelle, the 
orders were extracted using a direct sum method instead of the profile-weighted method used for 
the Ritter observations. Details of the reduction package can be found in Gordon & Mulliss (1997). 
The UT date and time, exposure length, and heliocentric Julian Date (HJD) for the observations 
are listed in columns 1-4 of Table |l[ 



3. Analysis 

From inspection of the spectra, we found that the lines of the two components of HD 197770 
were never fully separated. Figure [l] displays the C II A6578 A line for typical KPNO and Ritter 
spectra. The separation of the components in the two spectra are similar and the double-lined 
nature of this binary is clearly seen, especially in the KPNO spectrum. The lines are clearly 
resolved in the KPNO spectrum, but not in the Ritter spectrum. We identify the broad line 
as component A of this binary and the narrow line as component B. Since the lines were never 
fully separated, blending effects are always present. Therefore, we used TODCOR, the 2D cross 
correlation algorithm presented in Zucker & Mazeh (1994), to measure the radial velocities. This 
algorithm takes model spectra of the two components as input and finds the optimal radial 
velocities of the two components and their luminosity ratio. 

As results of TODCOR are sensitive to the assumed model spectra, we constrained the 
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Table 1. Observations and Radial Velocites 
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Fig. 1.— The C II A6578.03 A in the HJD 2448901.6 KPNO spectrum (bottom) and HJD 2450710.6 
Ritter spectrum (top). The resolutions of the spectra are 80,000 (bottom) and 25,000 (top). 
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model atmosphere parameters as follows. The effective temperature (T e ff) for HD 197770, derived 



using a reddening-free photometric index, is 21,000 ± 3,000 K (Gulati, Malagnini, & Morossi 



1989). We checked this T e g- by comparing the ultraviolet through V band dereddened spectral 



energy distribution (SED) of HD 197770 with PHOENIX LTE model SEDs (Aufdenberg et al 



1998). The SED of HD 197770 was taken from International Ultraviolet Explorer {IUE) and 



UBV (|Harmanec, Horn, fc Juza 1994 ) observations and dereddened assuming a Ry = 3.1 and 



E(B - V) = 0.58 flCardelli, Clayton, fc Mathis 1989|) . The unreddened SED was well fit by a 
model with a T e g = 21,000 K and log(g) = 4.0. In addition, since the eclipse depths are similar, 



the T eS of both stars must be close to 21,000 K QClayton 1996|) . We used a PHOENIX LTE model 



SED with a T cff = 21,000 K and a log (5) = 4.0 for both stars for the TODCOR algorithm. The 
vsini values were determined by iteratively running TODCOR and adjusting the v sinz values 
until the composite model spectrum matched, by inspection, the C II AA6578,6582 doublet in 
the HJD 2448901.6 spectrum (R ~ 80, 000). The strength of the model spectrum C II lines were 
reduced by 40% to match the observed strengths of the lines. This can be traced to the fact that 
the model spectrum had solar abundances and most B-type stars are underabundant compared 



to the Sun ( Snow fc Witt 1996 and references therein). The HJD 2448901.6 spectrum and the 
best fitting composite model spectrum is shown in Figure §. The best fit vsini values were 55 
and 15 km s _1 for components A and B, respectively. The TODCOR determined luminosity ratio 
(La/ Lb) was ~2. 

Using TODCOR, we were able to measure the radial velocity of the narrow, but not the 
broad, component in all of the spectra. We were able to accurately measure the broad component 
in only the 1992 October KPNO spectra. The 1992 June KPNO spectrum gave us radial velocities 
for both components, but neither was measured with great accuracy due to the low signal-to-noise 
ratio (S/N). The resolution and S/N of the Ritter spectra were not high enough to allow us to 
measure the velocity of the broad component (see Fig. |l|). As a result, the orbital motion of the 
narrow component was well sampled, but the orbital motion of the broad component was only 
measured 4 times, 3 times accurately. We fit the narrow component radial velocities using all 
but the HJD 2450360.6 measurement. This one radial velocity is quite deviant (4a from final 
fit, a = 2.8 km s -1 ) and a much better fit is achieved by rejecting this point. We fit the 3 good 
measurements of the broad component to determine Ka- For this fit, we assumed all the other 
orbital parameters were those determined from the narrow component fit, except wa = wb — 180°. 
The resulting orbital parameters are listed in Table [||. The values of m sin 3 i and a sin i were then 
computed from the orbital parameters. Figure || plots the measured radial velocities along with 
the radial velocities calculated from the orbital motion fits. Columns 5-9 of Table [l] give the phase, 
measured radial velocites (O), and observed minus calculated (O-C) radial velocities. 
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Fig. 2.— The C II AA6578,6582 doublet in the HJD 2448901.6 spectrum is plotted along with the 
best fit composite model spectrum. 
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Fig. 3. — The radial velocities for HD 197770 A & B are plotted phased to the period given in 
Table 2. The radial velocities calculated from the fits to the orbital motion are plotted as solid 
lines. 
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4. Discussion 



Since HD 197770 exhibits shallow eclipses (5V ~ 0.05), its inclination must be near 90° 
( Clayton 1996 ). Thus, the values quoted in Table [2] for msin 3 i are close to the actual masses 



of the components. Comparing an unpublished Pine Bluff Observatory (PBO) spectrum of 
HD 197770 to the spectra presented in Walborn & Fitzpatrick (1990), we find the spectral class 
to be Bl V-III or B2 III for both stars combined. Considering both the PBO spectrum and 
the T e ff value determined in §||, the most likely spectral type is B2 III. All lines in the spectra 
of HD 197770 are double, leading to the conclusion that both stars have similar spectral types. 
Assuming the radii of a B2 III star (12 Rq; Drilling fc Landolt 1997] ) for both components, the 



inclination of this eclipsing binary is > 73°. The masses of components A &; B are then < 3.3 and 
2.2 Mq, respectively. A normal B2 III star has a mass around 15 Mq ( prilling fc Landolt 1997] ). 



Thus, both components of this binary are undermassive for their given spectral types. This marks 
both components as evolved stars. 

HD 197770 lies in the Cygnus region on the edge of Cyg OB7 and Cep OB2. It seems to 
be very near and possibly on the edge of a molecular cloud/star formation region including 
Lynds 1036 and 1049. HD 197770 is associated with two IRAS point sources, IRAS 20420+5655 
and 20418+5700, but it is clearly non-stellar on the IRAS map at 60 [an. Figure |] displays the 
2° x 2° region centered on HD 197770 for the IRAS 12, 25, 60, 100 /im bands. The IR SED peaks 



at 60 /im implying warm dust in the immediate vicinity of the binary ( Gaustad fc Van Buren 



1993). In addition to the 60 /im shell with a radius of 14', the IRAS map shows the signature of 
an apparent bubble of cleared dust with a radius of approximately 0?6. This bubble is easiest to 
see in the 60 and 100 fim images, but is also present in the 25 //m image. The dark regions in 
the upper and lower right hand corners of the images mark the edge of a molecular cloud which 
has been mapped in CO ( pobashi et al. 1994| ). The molecular cloud wraps around three sides of 
HD 197770 roughly surrounding the evacuated area in the IRAS map. In particular, to the west of 
HD 197770 is a buried young stellar object (Lynds 1036) and to the east of HD 197770 is a pulsar 



( Dewey et al. 1985 ). So sequential star formation seems possible although the pulsar may be a 
background object. In any event, the HD 197770 binary appears to have formed on the edge of 
this cloud and cleared out an area around it. The association of HD 197770 with this cloud allows 



us to use the distance to the cloud (440 pc; Pobashi et al. 1994]) as an estimate of the distance to 
HD 197770. 

The existence of IR flux peaking at 60 fim at the position of HD 197770 and bubble of 
cleared dust surrounding HD 197770 imply that this binary has undergone at least two episodes 
of mass loss. A measurement of the luminosity of this binary would greatly help in determining 
its evolutionary stage, but the distance to the system is not known. HIPPARCOS gives a parallax 
of 0.52 ± 0.50 mas which results in a 3cr lower limit on the distance of 495 pc. This is consistent 
with the "guilt by association" distance given in the last paragraph. By assuming the distance 
to HD 197770 is 440 pc, Lb/ La = 2, and T e g = 21,000 K, we can estimate the luminosity and 
radii of the components of this binary. The resulting luminosities and radii are 9.3 x 10 3 L 



-9- 




20 40 60 80 

x [pixels] 

Fig. 4. — The IRAS 2° x 2° images centered on HD 197770 are displayed. The image intensity 
ranges are 1.2-2.5 MJy sr _1 for the 12 /im image (upper left), 4.2-4.95 MJy sr _1 for the 25 fim 
image (upper right), 4.0-6.0 MJy sr _1 for the 60 /im image (lower left), and 25-35 MJy sr _1 for the 
100 /im image (lower right). 
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(M v = -3.2) and 7.4 R for HD 197770 A and 4.7 x 10 3 L {My = -2.5) and 5.2 R for 
HD 197770 B. These radii result in an inclination for the binary of > 81° and, thus, masses of 
~3.0 and 2.0 M for HD 197770 A & B, respectively. 

We note the similarity between this system and the eclipsing binary v Sgr (m p sin 3 i = 2.5 M , 
m s sin 3 i = 4.0 M , & asini = 210 R ; Pudley &: Jeffery 1990 ) for which the primary is a 



hydrogen-deficient A-type supergiant and the secondary has a spectral type of B2 lb ( |Schoenberner 



fc Drilling 1983 ). According to Plavec (1973) and Schoenberner & Drilling (1983), the system 
is the result of type BB binary evolution, where the primary is now a post-AGB star. Uomoto 
(1986) has discussed the possibility that such systems are progenitors of lb-type supernovae. 

Two kinds of additional observations are needed of this binary. First, high resolution 
(R > 50,000) and high S/N optical spectra are needed spaced throughout the ~100 day orbit. 
These observations would confirm the K-amplitude of component A which is currently only based 
on three measurements. Second, imaging with an optical interferometer would give an accurate 
distance measurement when coupled to the results of this work. The Navy Prototype Optical 



Interferometer could perform this imaging when it is fully operational (Armstrong et al. 199 



Calculation of each star's luminosity would help in understanding their evolution stage. 



The Ritter observations were possible only with the help of the Ritter technician Bob 
Burmeister and the crack Ritter observing team. Team members contributing observations to this 
effort were Karl Gordon, David Knauth, Anatoly Miroshnichenko, Nancy Morrison, Christopher 
Mulliss, & Tracy Smith. Support for observational research at Ritter Observatory is provided 
by NSF grant AST-9024802 to B. W. Bopp and by The University of Toledo. This research has 
made use of the Simbad database, operated at CDS, Strasbourg, France. The IRAS data were 
obtained using the ISSA postage stamp server at IPAC. IPAC is funded by NASA as part of 
the IRAS extended mission under contract to JPL. K. D. Gordon & G. C. Clayton acknowledge 
support from NASA grant NAS5-3531. C. M. Anderson acknowledges support from NASA grant 
NAS5-26777. 
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Table 2. Orbital Parameters 



A B 


V [km s- 1 ] 


-15.31 ±0.32 


K [km s- 1 ] 


31.21 ±0.81 47.05 ±0.36 


e 


0.147 ±0.008 


w [°] 


69.76 ± 3.34 249.76 ± 3.34 


P [days] 


99.6918 ± 0.0223 


T G [days] 


41906.392 ± 1.975 


msin 3 i [Mq] 


2.89 ±0.08 1.92 ±0.09 


asinz [R©] 


60.8 ±1.6 91.7 ±0.7 



